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A simple and predictive model is put forward explaining the 
experimentally observed substantial shift of the glass transi- 
tion temperature, T g , of sufficiently thin polymer films. It 
focuses on the limit of small molecular weight, where geo- 
metrical 'finite size' effects on the chain conformation can be 
ruled out. The model is based on the idea that the polymer 
freezes due to memory effects in the viscoelastic eigenmodes 
of the film, which are affected by the proximity of the bound- 
aries. The elastic modulus of the polymer at the glass transi- 
tion turns out to be the only fitting parameter. Quantitative 
agreement is obtained with our experimental results on short 
chain polystyrene (Mw = 2 kg/mol), as well as with earlier re- 
sults obtained with larger molecules. Furthermore, the model 
naturally accounts for the weak dependence of the shift of T g 
upon the molecular weight. It furthermore explains why sup- 
ported films must be thinner than free standing ones to yield 
the same shift, and why the latter depends upon the chemi- 
cal properties of the substrate. Generalizations for arbitrary 
experimental geometries are straightforward. 



INTRODUCTION 

To explain the experimental observation that thin 
polymer films melt at temperatures well below the glass 
transition temperature of the bulk polymer, T° is 
one of the major current challenges in the theory of ho- 
mopolymers. At first glance, one might expect such a 
behavior, due to the impact of the finite size geometry of 
a thin film upon objects as large as polymer molecules. 
In fact, for molecular weights larger than about M\y = 
300 kg/mol, geometry effects have been recently shown 
to play a dominant role, resulting in a linear depen- 

For smaller 



dence of T g upon the film thickness, h 
molecules, however, the problem is conceptually more in- 
tricate, since a noticeable reduction of T g can be observed 
at film thicknesses orders of magnitude larger than the 
radius of gyration of the molecules. Furthermore, the 
shift Tg — T g (h) becomes strongly nonlinear and largely 
independent of molecular weight jj|7), indicating that 
this regime must be governed by a different mechanism. 
A fundamental understanding of this effect would be of 
great interest not only for polymer physics, but poten- 
tially also elucidate the physics of the glass transition in 
a larger class of systems, since it comes into play at small 
molecular weight. 



Many attempts have been undertaken to explain these 
observations, mostly by considering microscopic models 
of the inner structure of the films. D. Long and F. 
Lequeux have envisaged the freezing of the film as a per- 
colation of rigid domains mediated by thermal fluc- 
tuations H . Other models assume a layered structure of 
the film, with a particularly mobile region close to the 
free surface of the film pL Within this framework, J. 
Forrest and J. Mattsson [7fl have recently been able to 
achieve quite impressive accordance with the experimen- 
tal data [§[0]. Their model makes use of the so-called 
cooperativity length, £(T), which plays a mayor role in a 
whole class of theoretical concepts of the glass transition. 
The only drawback is that there is yet no well established 
theory of £(T). Furthermore, as a consequence of the 
two-layer structure of the film inherent in the model, it 
is not completely clear why there should not be two glass 
transitions, rather than a single one shifted in tempera- 
ture. 

As a complementary approach, computer simulations 
of polymer films with free surfaces have recently been car- 
ried out, and most of the experimental findings were re- 
produced qualitatively ]l0| . However, the polymer chains 
in these simulations were only 16 monomer units long, 
very much shorter than those used in the experiments. 
Before attempting to explain the effect theoretically, it 
is therefore worthwile to study experimentally the be- 
havior of polymers with short chain length, both to ease 
comparison with simulation and to explore the range of 
validity of the apparent independence of the shift of T g 
on molecular weight. 

Before going into the details of our study, let us take a 
break and summarize the main questions to be answered. 

1. What is the principal mechanism responsible for 
the reduction of T„ in thin films of low molecular 
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weight {Mw <300kg/mol) polymers 



2. Down to how small molecular weight is this mech- 
anism valid? 

3. Why is the effect stronger in free standing films 
than in supported ones [gj? 

4. Why is there no significant dependence of this effect 
on molecular weig ht @? 

5. Why does the effect depend upon the chemical 
composition of the substrate for supported films 

Ml? 
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6. Why is there sometimes an increase of T g in thin 
films, instead of a reduction [p]-p^[? 

We will try in the present paper to give answers to these 
questions, or at least show in which direction answers 
might be found, on the basis of a novel, quite simple 
model which rests mainly on the viscoelastic eigenmodes 
of the films. We restrict the discussion to the case of 
polystyrene (PS), since this is the most thoroughly stud- 
ied polymer in this context. Furthermore, it is particu- 
larly well suited for comparison with theoretical models, 
since (atactic) polystyrene does not show any propensity 
to crystallization. In other polymers, which might crys- 
tallize at least in part of the film, modelling would be 
exceedingly difficult, and the main mechanism could well 
be obscured in experiments. 

EXPERIMENT 

Let us first explore the range of validity of the afore- 
mentioned effect of reduction of T g , as to the molecular 
weight of the polymer. We have investigated the glass 
transition in thin supported films of PS with a molecular 
weight as small as 2 kg/mol (« 20 monomer units). It 
was purchased from Polymer Labs (UK) with a polydis- 
persity index M\v/Mn = 1.05, the radius of gyration is 
1.3 nm. Effects from the molecular geometry are thus 
expected only for films of few nanometers thickness. The 
films were spin cast from toluene solution onto silicon 
wafers (Silchem GmbH, Freiberg/Germany), which were 
previously cleaned by ultrasonication in acetone, ethanol, 
and toluene, subsequently. Residual organics were re- 
moved with a 1:1 mixture of H2SO4 with H2O2, and the 
substrates were thoroughly rinsed with hot millipore wa- 
ter afterwards. Films were investigated with thicknesses 
ranging from 4 to 160 nm. The roughness of the free 
surface of the films was less than 0.2 nm, as revealed by 
scanning force microscopy (SFM). 

The glass transition temperature was determined in 
two different ways, depending on film thickness. The 
standard procedure of monitoring the thermal expansion 
of the film via ellipsometry, as introduced by Keddy and 
Jones J!]], was used for film thicknesses down to 9.6 nm. 
Fig. |l|a shows a typical run. The data were reproducible, 
irrespective of being taken during heating or cooling. 
Typical heating or cooling rates were 2 K/min. We plot- 
ted the changes in refractive index and thickness of the 
film in fig. [l|b. The solid line represents the Clausius 
Mosotti relation. Obviously, there is good agreement, 
suggesting the absence of any loss or degradation of ma- 
terial. This was found invariably for all samples. 

In fig. H we plotted the thermal expansion coefficients 
found above and below T g (h), which coresponds to the 
kink in fig. [lja. At large film thickness, T g obtained in this 
way approached 327 ± 1 K, which is consistent with the 



temperature at which macroscopic melting is observed in 
the bulk (T°) for PS with this chain length. The shaded 
areas represent published data of the expansion coeffi- 
cients JT3]. Good agreement is found both above and 
below T g . Small deviations at very small film thickness, 
as might be concluded from our data, have been reported 
before 0. 

For films thinner than 9 nm, we determined the melting 
behaviour of the films by observing the buildup of am- 
plified thermal fluctuations (spinodal dewetting) p^| , ^6[ . 
These processes were monitored by SFM with in-situ 
heating. In order to speed up the experiments to a fea- 
sible time scale, dewetting was observed at temperatures 
close to Tg. As an example, we show in fig. || the tem- 
poral evolution of the Fourier transform (spatial power 
spectrum) of the surface topography. The hallmark fea- 
ture is the clear peak representing the fastest growing 
mode. The inset shows the peak intensity as a func- 
tion of time on logarithmic scale. Exponential growth is 
clearly observed up to a rather well defined time at which 
coalescence of holes sets in. 

From the slope of the straight line in the inset, and the 
known effective interface potential jl6| , the viscosity r\ of 
the polymer can be determined. Measuring the viscosity 
as a function of temperature, we found that this obeyed 
a Vogel-Fulcher law, with the Vogel-Fulcher temperature 
shifted by a certain amount ATvf, which depended on 
the film thickness. We identified A7V^(/i) with the shift 
in glass temperature. In this way, the apparent glass 
transition temperature of particularly thin films was in- 
ferred from the spinodal dewetting experiments by set- 
tingT g (h)=T°-AT VF (h). 

Our experimental results are shown in fig. ^ as the 
full symbols. The circles represent the thermal expan- 
sion measurements, the squares were obtained from the 
spinodal dewetting experiments. As one can clearly see, 
the glass transition temperature is substantially reduced 
for all films thinner than about 50 nm. The solid line 
represents the function 

Tg = T°(l + ho/h)- 1 (1) 

This form has been shown before to account well for the 
data obtained by others for larger molecular weight films, 
if h — 0.68 nm was assumed for PS p7| . Within ex- 
perimental scattering, our data exhibit indeed the same 
dependence of T g (h) in the full range of film thickness 
explored. It is remarkable that the data from both the 
thermal expansion and spinodal dewetting measurments 
are well fitted by the same curve. This tentatively corrob- 
orates the procedure of obtaining T g for very thin films 
(squares) as discussed above. 

Although our polymer chains are roughly by a factor 
of 50 shorter than those investigated before, we obtain 
ho = 0.82 from the fit, which is quite close to the above 
value. This confirms the weak dependence (if there is a 
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significant one at all) of the reduction of the glass transi- 
tion temperature on the molecular weight of the polymer, 
down to a molecular weight as small as 2 kg/mol. 

We can thus state that the reduction of T„ in thin films, 
as described phenomenologically by eq. (pi), is observed 
in a huge range of molecular weight, from molecules as 
small as 2 kg/mol up to a few hundred kg/mol. This is 
in accordance with molecular dynamics simulations [jlO| , 
and rules out mere finite size effects on the individual 
coils as the main cause of the reduction of T g in this 
regime. Note that the radius of gyration of our poly- 
mer is only 1.3 nm, while T g is significantly reduced at 
a film thickness of 50 nm already. The fitting parameter 
used in eq. ([!]), h , changes only by about 20 % in this 
range, confirming that the dominant mechanism which is 
responsible for this effect cannot depend strongly on the 
molecular weight of the polymer. 



EIGENMODE SPECTRUM OF THE FILMS 

Since it is clear that the behaviour displayed in fig. || 
can in no way be attributed to the geometrical impact of 
the finite film thickness upon the microscopic conforma- 
tion of the individual chains, we present here an approach 
to the problem which intentionally makes as few refer- 
ence as possible to the molecular structure of the film. 
The latter is accounted for merely by the strain in the 
polymer, i.e., the deviation of the local average gyration 
ellipsoids of the molecules from a sphere. By 1 local! we 
mean a volume much larger than the volume of the back- 
bone of a single molecule, but with a lateral dimension 
much smaller than the film thickness. 

Such a deviation from the equilibrium conformation 
(i.e., from a Gaussian coil, if self-avoiding is neglected), 
which may be viewed as an entropy fluctuation, can de- 
cay either by self diffusion of the individual molecules, 
or by some center-of-mass rearrangement (i.e., flow) of 
the melt. It is clear that close to a surface, such rear- 
rangements are much easier to accomplish due to cou- 
pling of the capillary waves on the free surface to the 
bulk flow of the polymer. This coupling is effective down 
to a depth comparable to the wavelength of the capillary 
modes, which may be much larger than the coil size of 
the molecules. Thus we are provided with a mechanism 
which affects the motion of the polymer molecules, and 
naturally can act over distances which are large as com- 
pared to molecular dimensions. In order to explore the 
possible relevance of this mechanism for the reduction 
of T g in thin films, we have to consider the eigenmode 
spectrum of the (viscoelastic) polymer film. This can be 
discussed with all possible boundary conditions at the 
substrate, such that films with strong slip along the sub- 
strate, grafted films, or free standing films, may as well 
be treated within the same framework. 



The spectrum of a viscoelastic thin film can be ob- 
tained in a straightforward manner by combining stan- 
dard theory of elasticity Jl8| and hydrodynamics (l9j in 
the limit of small Reynolds number (Stokes dynamics). 
The equation of motion reads 



7] T] 
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where E is Young's modulus, r\ is the viscosity, and u>o 
is the Rouse rate of relaxation of the individual chains 
into their equilibrium configuration ]20| , pl"| ] . <fi is a vector 
field related to the strain tensor, S. For the sake of clar- 
ity, we restrict our discussion to a simple model, taking 
into account only the single intrinsic relaxation rate loq, 
as opposed to more general treatments p2j . Deviations 
of the real polymer from this simple behaviour will be 
introduced a posteriori farther below. 

If we restrict the discussion to one lateral (x) and one 
normal (z) coordinate, (f> — (^x^z) is defined via 



S = 



\{dx<t>z + d z (j) x ) 



\{d x (j) z + d z (f>x) 



(3) 



Finally, p is the pressure field. 

For harmonic excursions of the free surface, ((x) = 
Coexpjiga; — tot}, eq.(|^) has solutions 

(f) x = [1 + (h + q~ 1 )a(q)] coshqz — q~ 1 a(q) sinhgz 

4> z — [1 + ha(q)] smhqz — za(q) coshgz (4) 

where for the function a(q), we find 



a{q) = I2J e^-l + qh 



(5) 



for free standing films (symmetric modes) as well as for 
supported films with full slippage (zero friction). For 
supported films with some friction at the substrate, the 
expressions are of similar form, but considerably more 
complicated, and will not be discussed here. At the free 
surface, we used the standard boundary condition of zero 
tangential stress, andp = —ad xx (, where a is the surface 
tension of the polymer. Note that for free standing films, 
h is defined as half the film thickness. 

For the relaxation rates of the modes, we get 
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Since it is only the 'fast' modes which contribute appre- 
ciably to the reduction of the glass transition tempera- 
ture (see below), we consider only the upper branch of 
eq. (||). Observing that luq << —, this is given by 
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as a very good approximation. The dependence of the 
wave number, q, stems solely form the coupling to the 
capillary waves on the free film surface. The first major 
assumption of our model is that T g (h) is determined by 
the spectrum of the viscoelastic eigenmodes of the film 
as given by eq. (0). 



THE FREEZING MECHANISM: MEMORY 
EFFECTS 

The second major assumption is that the physical 
cause for the melting or freezing of the film, respectively, 
are memory effects in the polymer material. These are of 
course not included in the linear theory discussed above, 
and may be formulated in a generic way by means of a 
suitable memory kernel, as used in a class of theoreti- 
cal models of the glass transition, called mode coupling 
models |p3| , p4j| . Within this framework, memory effects 
are taken into account by inserting a convolution integral 
with the memory kernel m{(f>(t)} — ai(/>+a20 2 + &3</> 3 + ... 
in the otherwise linear differential equation of motion of 
the modes considered p5[ . In these models, cj> usually de- 
scribes density fluctuations. In contrast, we consider the 
material to be essentially incompressible, and (f> denotes 
here the strain in the polymer material, as defined above. 
It thus describes the local state of the material, similar 
to what the density does for simple glass forming liquids 

His 

In our case, the equation of motion (2) is of first order 
in time, hence we have 



cj)' + uj{q)cj) + / m{(f>{T)}(f>' '(t - r)dr 







(8) 



as the mode coupling equation. This type of equations 
has been thoroughly analysed [^,^,^7) in relation to 
the microscopic physics of the glass transition (to which 
we do not refer here), as well as to large scale degrees 
of freedom [^8|. For density fluctuations in glass form- 
ing simple liquids, it was found that the coefficients a, 
of the memory kernel vary concurrently with tempera- 
ture, and that upon crossing a certain border in the space 
spanned by the c^, the system freezes into a nonergodic 
state |m| . The existence of such a freezing transition has 
been found to be largely independent of the precise form 
of the memory kernel. In fact, most of the features of a 
glass transition may be well represented in what has be- 
come known as schematic models, which are simple mode 
coupling equations not referring in their memory kernels 
to the microscopic physics of the system under study. 

It is now worthwhile to contemplate on possible mem- 
ory effects in a polymer melt, as to their scaling with 
temperature. Let us first consider a polymer molecule 



in equilibrium, forming a more or less Gaussian coil. If 
this is elongated by straining the polymer melt, to what 
extent will it memorize this process after the strain is re- 
leased? In the strained state, the molecule will relax to 
some extent. However, this relaxation will not proceed 
homogeneously along the molecule, since the activation 
energy, U, for local rearrangements will depend upon the 
local topological environment. 

In order to discuss the relaxation behaviour of the lo- 
cal molecular geometry, let us define a local geometric 
exponent of the coil, //, by the relation 



< A | r |>= (As) 



(9) 



where < A | r | > is the typical distance travelled in space 
upon moving alongside the polymer chain by a distance 
As. In completely stretched parts of the chain, /.i = 1, 
whereas in equilibrated domains, p. = [lequUib. ~ 0.5. 
The variation AU of U(s) determines the variation of the 
relaxation rate of \i towards fi e quiUb. ■ When the strain is 
released, remnants of this variation will remain, and thus 
represent a memory of the strain. This holds as long 
as the strain/release process is fast as compared to the 
equilibration time, w^ 1 ■ For the modes to be considered 
here, this is well fulfilled. When AU << kT, as can be 
safely assumed, these remnant variations in p(s) scale 
as AU/kT, thus we can conclude that memory effects in 
the polymer scale as 1/T. As a direct consequence, we 
can replace the memory kernel m{4>(t)} of eq. (||), the 
coefficients ai of which are temperature dependent, with 
M{4>}/T, where M is now independent of temperature. 

By proper normalizaton of time, eq. (^|) can now be 
rewritten as 



M{(j)(T)}(f>'(t-T)dT =0 (10) 



Tu(q) 



From this equation, we see directly that the modes with 
the largest relaxation rate freeze at the lowest tempera- 
ture. Furthermore, the precise form of M, which might 
be derived from a detailed analysis of the non-equilibrium 
dynamics of the polymer molecules, does not need to be 
considered here any further. It is only required that it 
belongs to the class of kernels which yield a freezing tran- 
sition for <f> at all- 
Let us now turn back to the eigenmodes of the 
film. The relaxation rates, as given by eq. (^), are 
monotonously increasing with q. However, modes with 
q much larger than the inverse film thickness, h , do 
not penetrate appreciably into the film, such that only a 
small fraction of the material takes part in these modes. 
Hence we are led to considering chiefly the modes with 
q « h , since these are the highest frequency modes 
comprising all of the film material. 



Setting thus T^ft,- 1 ) 



const, as suggested by 



eq. (|10|), we directly arrive at a simple formula for the 
glass transition temperature: 
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T g (D) = T°ll 



1.16 a 



h(E + rjujo) 



(11) 



This is precisely the form of eq. |l]), and fig. || shows that 
it describes our data very well. 



DISCUSSION 

The quantity rju>o is on the order of a few kPa and can 
in general be safely neglected against E. The characteris- 
tic length scale which appears here, and which was called 
ho in eq. ([j]), is thus the ratio a/E, with the surface ten- 
sion a = 31 mN/m for PS. It follows that the elastic mod- 
ulus determining the dynamics of the relevant modes, E, 
is the only physical fitting parameter in the model. From 
h{) = 0.82 nm, as obtained from the fit of our model, we 
find for the elastic modulus E « 44 MPa, which is, on 
logarithmic scale, right in between the modulus of the 
frozen material (a few GPa) and the modulus just above 
T g (about 300 kPa) . On the basis of the model discussed 
above, let us now try to formulate answers to the ques- 
tions asked in the introduction, one after the other. 

1. What is the main mechanism responsible for the 
reduction of T g in thin films of low molecular weight 
{M-w <300kg/mol) polymers? The physical picture 
which emerges from our model is that as the tempera- 
ture is increased, melting proceeds as the fastest mode 
involving all of the film material escapes from its frozen 
state and fluctuates. This proceeds with the help of the 
capillary waves on the free film surface, which ease en- 
tropy fluctuations in the polymer melt. The concomi- 
tantly increased motion of polymer chains reduces the 
effective viscosity also for modes with smaller q, which, 
as a consequence, are sped up (77 appears in the denom- 
inator in eq. (0)!) and melt in turn. In this way, the 
film finally melts at all length scales. It is illustrative to 
note that the softness of the spinodal modes, as used for 
the determination of T g at small film thickness, does not 
contradict our model: these are modes with q << l/h, 
and are thus not the fastest modes, as considered in the 
melting mechanism. 

2. Down to how small molecular weight is this mecha- 
nism valid? We found that in agreement with simulation 
results p0[ , the effect is present even for molecules as 
short as about 20 monomer units. This is well below 
the limit of entanglement. It thus spans a wide range 
in molecular weight, over more than two orders of mag- 
nitude. It is precluded, or obscured, by other effects at 
molecular weights in excess of about 300 kg/mol ||. 

3. Why is there no significant dependence of the ef- 
fect on molecular weight? The physical property which 
solely determines the thin film behaviour according to our 
model is the elastic modulus, E. As it is well known, E 
displays no marked dependence on the molecular weight 
P4E29II , such that on the basis of our model, T g (h) is 



expected as well to be largely independent of molecular 
weight. To be precise, E is slightly less for smaller molec- 
ular weight than for larger fll|,^J , such that according to 
eq. (pi]), ho should be larger for smaller molecular weight. 
In fact, we obtained ho = 0.82 nm for our very short 
molecules, whereas for larger molecules, ho = 0.68 nm 



was found |17|. The observation that the relevant length 
scale, hg, decreases with increasing molecular weight is 
another evidence against geometrical effects on individ- 
ual chains to be relevant in the regime discussed here. 
The exact physical significance of the somewhat arbitrary 
'choice' of E, which may be viewed as E at T g , is to be 
investigated in further studies. 

4. Why is the effect stronger in free standing films 
than in supported ones? We have mentioned that the 
eigenmodes of supported films of thickness h with no fric- 
tion are identical to those (the symmetric ones) of free 
standing films of thickness 2h. Thus our model predicts 
that the effect in a supported film is just as large as in a 
free standing film of twice the thickness. This is indeed 
in accordance with experimental observation || . 

5. Why does the effect depend upon the chemical com- 
position of the substrate for supported films? It is clear 
that the amount of friction of the film material at the sub- 
strate changes the eigenmode spectrum of the film, such 
that a dependence on the chemical composition of the 
substrate is indeed expected on the basis of our model. 
For a quantitative comparison with experimental data, 
however, the friction coefficients would have to be de- 
termined for the systems investigated. It should be well 
noted that this also induces some uncertainty in the nu- 
merical value of E derived from our fit (cf. fig. ^), since 
we have not characterized the friction of the films on the 
substrate. 

6. Why is there sometimes an increase of T g in thin 



films ]11| |13J? We can try to give a rough idea to what 
extent effects like this may be accounted for within our 
model, by including the interaction of the free surface 
with the substrate via long range forces. This is described 
most conveniently by replacing aq 2 by aq 2 + , where 
V(h) is the effective interface potential of the film due to 
long range forces ]30|j3l]| . For unretarded van der Waals 
forces, V(h) = — AjYlTth 2 , where A is the Hamaker con- 
stant. If A > 0, is negative, such that the van der 
Waals forces tend to destabilize film, giving rise to spin- 
odal dewetting for sufficiently thin films. In any case, 
a positive A will reduce the relaxation rate of all of the 
modes, including the one at q = l/h, and thereby tend 
to increase the glass transition temperature. Quantita- 
tively, this effect is obtained by accordingly replacing a 
in eq. ( pi] ) by a — 2^2- As it turns out, T g attains a 
minimum at h m i n = 3^ A/2ira, and increases sharply 
for smaller thickness. For A — 2.2 x 10~ 20 J, represent- 
ing PS on silicon oxide p2]| , this is at /i mill w 1 nm. An 
increase of T g above T? at a film thickness significantly 



5 



larger than 1 nm can be explained only with an unphys- 
ically large Hamaker constant. 

However, one should anyway be cautious with trying 
to interpret an increase of T g in thin films in the frame- 
work of our model if working with polymers other than 
PS [jlO] [l3|]. It is not at all clear what impact a possible 
crystallization (which is not present in atactic PS) can 
have on T g (h). Consequently, possible model systems for 
studying the mechnaism discussed in the present paper 
have to be carefully chosen to ensure the absence of any 
pronounced inner texture of the film. Our model can 
be appropriate only for 'structureless' films, and might 
break down for films with a pronounced layered texture 
like, e.g., a Langmuir-Blodgett film [j§3) , or for other poly- 
mers which may crystallize at least in some part of the 
film. 



OUTLOOK 

It is tempting to apply the view developed here also 
to the surface of a bulk polymer sample, for which the 
spectrum is obtained setting a = q/2. As it is readily 
seen, one should expect surface melting at the polymer 
surface, down to a thickness of h sm = h Tg(Tg — T) _1 . 
This is very much along the lines of first ideas which had 
been put forward to explain the observed reduction of 
T g in thin films [EJ. To the best of our knowledge, there 
is yet no conclusive evidence in favor or in disfavor of 
polymer surface melting. 

To confirm the theoretical concept put forward here, 
an ab-initio derivation of the memory kernel, m{0}, from 
the statistics and microscopic transport properties of 
the individual molecule would be desirable. Further- 
more, it is important to elaborate on the significance of 
Young's elastic modulus, E, at the glass transition. This 
was determined experimentally using our model to be 
E = 44 MPa, but it is not yet clear what physical prin- 
ciple distinguishes this value. A complete theory of the 
glass transition on the basis of the proposed model would 
be rather cumbersome, since it involves two strongly non- 
linear mechanisms mutually affecting each other: freezing 
by the mode coupling mechanism, and non-Newtonian 
effects on the viscosity. Our model furthermore sug- 
gests that the glass transition in a homopolymer may 
be viewed as a mode-coupling-model freezing of its vis- 
coelastic bulk modes. This needs not to contradict more 
classical views of the glass transition in homopolymers, 
but might serve as an interesting, and potentially useful, 
alternative approach. 
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FIG. 1. a) Ellipsometric dilation measurement of a thin 
PS film. Typical heating/cooling rates were 2 K/min. b) The 
change in film thickness vs. the change in the optical dielec- 
tric constant, as measured by ellipsometry for various temper- 
atures during a heating/cooling cycle. The good agreement 
obtained with the Clausius-Mosotti relation represented by 
the solid line demonstrates the absence of significant loss or 
degradation of material. 

FIG. 2. Thermal expansion coefficients obtained above and 
below T g . The shaded areas represent values taken from the 
literature Q. 

FIG. 3. Temporal evolution of the power spectrum of spin- 
odal dewetting undulations on a PS film exhibiting unstable 
capillary waves at its free surface. Curves of equidistant times 
are superimposed to show the gradual increase of the main 
peak. The inset shows the exponential growth of the latter, 
from which the viscosity of the film can be inferred. 

FIG. 4. The glass transition of thin films of 2 kg/mol 
polystyrene, as determined from thermal expansion (circles) 
and from the growth of spinodal waves (squares). The solid 
curve represents our model, which has the elastic modulus 
governing the dominant modes as its only fitting parame- 
ter. It furthermore corresponds to what was found for larger 
molecular weight PS before jfjj|. Top: linear scale. Bottom: 
logarithmic scale, showing more details at small film thick- 
ness. 
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